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Abstract 



m 
o 

^ I We compare the speed with which a sexual, respectively an asexual, population is 

able to respond to a biased selective pressure. Our model focuses on the Weismann 
, hypothesis that the extra variation caused by crossing-over and recombination during 

^ I sexual reproduction allows a sexual population to adapt faster. We find, however, that 

CN ■ the extra variation amongst the progeny produced during sexual reproduction for most 

model parameters is unable to overcome the effect that parents with a high individual 
fitness in general must mate with individuals of lower individual fitness resulting in a 
■ moderate reproductive fitness for the pair. 

o 

Q"^ , Keywords: mode of reproduction, cost of sex, speed of adaptation. 
O . 

2 . Short title: Sex and reproductive fitness. 

g : 1 Introduction 

Q I It is generally assumed that an important contribution to overcoming the evolutionary cost 
CJ ' 

of sexual reproduction is the greater genetic variation among the progeny of sexual parents, 
^ I caused by crossing-over and recombination, compared to asexually reproduced offspring. 

The greater variation in the sexual population is assumed to produce more choice for natural 
selection and might thereby enable swifter adaptation. This is the hypothesis of Weismann 
|Burt, 2000| . Several recent experimental studies have been interpreted as supporting Weis- 
mann's hypothesis [Colegrave, 20021 [Colegrave et al., 20021 |Kaltz and Bell, 2002| . These 



experiments observe a larger variation in the reproductive fitness (rate of reproduction) 
when sexual reproduction is involved and this produces a faster adaptation of the sexual 
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population. Thus it might seem that Weismann's suggestion has verified and fully under- 
stood. 

However, we believe that there are still issues to be clarified concerning the details 
of the mechanisms behind Weismann's hypothesis. The question we discuss here is the 
following: It is evident that crossing-over and recombination may lead to greater genotypical 
variation amongst the progeny but what matters for selection is the differential excess fitness 
(reproduction rate) of sexual reproduction compared with asexual reproduction. How is the 
greater genotypical variation amongst individuals able to lead to a larger reproduction rate 
of the pair of parents'? The essential point is that in a sexual population the individual of 
highest fitness typically will have to mate with an individual of lower fitness. 

We use here a simple schematic mathematical model of evolution in genotype space 
and study the speed with which the population is able to climb a fitness gradient. Sex- 
ual individuals can mate at random with any other individual of the same species, which 
we, according Mallet [Mallet, 1995 genoptypic cluster definition, interpret as individuals of 



sufficient similar genotype. Although our model is simple, it is nevertheless representative, 
as far as testing Weismann's hypothesis, of any approach in which the fitness of the sex- 
ually reproducing pair is obtained by combining measures of fitness of the two individuals 
involved in the reproductive event. We find that the sexual population generically exhibit 
a greater spread in genotype space compared with the asexual population. However, only 
in extreme cases does this greater variation in genotype space lead to a larger variation in 
the reproduction rate and thereby to a superior adaptability for the sexual population. 

2 The Model 

Genotype space: We consider a population of individuals each characterised by their specific 
genotype determined by a position vector S in genotype space S. For simplicity we model 
5 by a two dimensional x Ly grid, i.e. S = {Sx,Sy), where 5*^. and 5*^ are integers: 
Sx = 1,2, . . . , and Sy = 1,2,..., Ly, with periodic boundary conditions in the y-direction. 
We have used = Ly = 100 for the study reported here. The number of individuals at 



position S at time t is given by n{S,t). 

Reproduction rate: The genotype determines directly the reproduction rate, or the 
reproductive fitness, 0(S,t) in the following way 

The statistical weight W(S,t) is given by 

W{S,t)=aS,- fiN{t). (2) 

Here the first term represents a linear increase [a > 0) in the reproduction rate along the 
x-direction in genotype space. In the second term N{t) = X]s^(S,^) represents the total 
number of individuals at time t, and /i determines the carrying capacity of the system. For 
fixed N(t) the fitness landscape corresponding to W{S,t) in Eq. Q is a smooth upwards 
tilted plane. Assume for a moment that the entire population is placed on one location in 
genotype space, say S", and assume further more that adaptation cannot take place (zero 
mutation rate). A stationary population size is then determined by the balance between 
the reproduction [A^ —>■ N + N(f){S)] and ensuing annihilation of individuals [(1 + (/)(S))A^ — *• 
(1 + 0(S))A^(1- pkiu)], leading to 0(S) = pkui/ (1 - Pkui)- From this expression we find the 
population size, for a given carrying capacity /i, by use of Eq. ([T)) and ©, 

When mutations are allowed the population can adapt (done by increasing the component 
of the position in genotype space) to the environment and thereby increase the population 
size for a given value of the carrying capacity parameter fi. We have used the following 
parameter values: a = In 3/45 and fi = In 3/450 and p^m = 0.2 which leads to the linear 
dependence (A^) = 10{Sx) +450 of the average number of individuals (A^) on the population 
averaged x-coordinate (Sx) in genotype space. 

In order to be able to represent a more rough landscape we choose a fraction of 
genotype positions Sh at random and denote them fitness holes. The y-direction in the 
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genotype space is essential in this case as this extra dimension allows a population to bypass 
a fitness hole. These sites are infinitely unfit corresponding to (f){Sh) = 0, i.e. individuals 
with a genotypical composition equal to one of the Sh locations cannot produce offspring. 
The fitness landscape produced by this procedure resembles Gavrilets' notion of the holey 



Time step : The system is simultaneously updated. First all individuals are allowed 
to reproduce by the relevant probability to be described below. After reproduction annihi- 
lation occurs, during which individuals are removed from the system with probability Pkui- 
This probability is constant in time and equal for all genotypes. 

Asexual reproduction: Let us now describe how reproduction is represented in the 
model. We begin with the asexual case. All N{t) individuals are allowed to produce one 
offspring with the genotype specific probability po//(S) = 0(S, t). The offspring will be of the 
same genotype S except when mutations occur, which happens with a constant probability 
Pmut, to be thought of as the genomic mutation rate per reproduction event. The mutant is 
placed at random on one of the four nearest neighbour sites. 

Sexual reproduction: An individual can mate with any other individual if the two are 
sufficiently similar to belong to the same species. We take this to mean that two matting 
individuals must share a certain degree of genotypical similarity, cf. Mallet's genotypic 
cluster species definition Mallet [Mallet, 1995| . We consider two individuals to be of the 
same species if they belong to genotypes which are separated by less than a distance dmax 
in genotype space. 

First a mate of genotype Smate is found for each individual i of genotype Sj for 
2 = 1,2,..., N{t), This is done by randomly choosing the mate amongst all the individuals 
present on genotypical positions less than a distance d^ax away from Sj. In this way we 
ensure that only individuals belonging to the same species mate and that there is no pref- 
erential mating within the species. Each of the N{t) pairs will produce offspring with the 
respective probability 



fitness landscape [Garvilets, 19991 [Gavrilets et al., 1998| [Gavrilets et al., 2000 . 




(4) 
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The variation amongst the offspring induced by crossing-over and recombination is repre- 
sented by placing the offspring on a site So// selected randomly with uniform probability 
within an ellipse with eccentricity e and major half axis of length a = [l/2+|Si — S2| + l/2]/2e 
placed symmetrically around Si and S2, see Fig. ^ The likelihood with which the offspring 
are allowed to end up outside the region in genotype space between the two parents is charac- 
terised by what we call the recombination factor r = 1/e. The larger the value of r the more 
likely it is to find the offspring to the right (or to the left) of both parents, see Fig. ^ In 
the special case where both parents are of identical genotype, i.e. Sj = Smate, the offspring 
is given the same genotype as the parents, i.e. So// = Sj. 

Mutations can occur with the same probability Pmut as in the asexual case. When a 
mutation occur the offspring is moved to one of the nearest neighbour sites of the position 
So//- 

We emphasise that the specific function in Eq. is inessential. We obtain qualita- 
tive identical results when we instead of Eq. Q use other functional forms for po// which sat- 
isfies the seemingly reasonable requirement that Poff is a number between the individual fit- 
nesses 0(Sj) and (f){Smate)- for example one can replace Eq. (jlj) bypo// = [0(Si)+0(Smate]/2. 

It is also important to point out that although the population grows monotonically 
in size as it moves to the right in genotype space, the average reproduction probability 
(po//(S)) remains constant. The reason for this is that the increase in fitness caused by the 
increase in aS^ in Eq. Q is compensated by the carrying capacity term —fiN{t) in this 
equation. 

3 Results 

In order to compare the efficiency with which the two different types of reproduction respond 
to a selective pressure, we compare the evolution under the same conditions (i.e. same a and 
/i in Eq. (0)) of a solely sexually reproducing population with a solely asexually reproducing 
population. In both the asexual and sexual case five different genotypes are initiated with 
100 individuals each. They are chosen so that the average S'a,.-coordinate (Sx) = 5. The five 
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positions are S = (4, Ly/2), {5,Ly/2), {6,Ly/2), {5,Ly/2 + 1), {5,Ly/2 — 1). This gives an 
initial reproduction probability for all the 500 individuals of about Pojf = 0.252. The exact 
number of starting genotypes does not affect the nature of the movement through genotype 
space. 

We monitor the time it takes the population's centre of mass (S) to move in genotype 
space from (S) = 5 to (S) = 85. During this motion we find that (S) to very a good 
approximation increases linearly with time. 

Under the evolutionary dynamics the population moves to the right. Although the 
dynamics is diffusive in nature the population remains confined to a rather narrow region 
along the Sx and Sy axes while the population gradually moves to lager Sx values. This 
is shown as the peaks in n{S,t) in Fig. |21 The dispersive action of the diffusion in geno- 
type space produced by the mutations, and the crossing-over effect in the sexual case, is 
counteracted by the effect of —^N(t) in Eq. The reproduction probability poff{S) of 
individuals left behind, at genotype positions with small Sx values, will decrease with time 
because the term —^N{t) increases as a result the reproductive activity of individuals with 
higher Sx values. As soon as the reproductive probability of a site S becomes permanently 
smaller than the killing probability, Poff{S,t) < Pkuu the population on that site is bound 
to go extinct. 

We now summarise the behaviour of the model. First we consider the case without 
any fitness holes in genotype space. The sexually reproducing population is found only 
to adapt more rapidly than the asexually reproducing population for substantial variation 
amongst the progeny or low mutation rate. This is illustrated in Fig. IHl For a given mutation 
rate, the data points indicate the threshold value, rthri of the recombination factor where, 
for increasing r, sexual reproduction becomes more effective than asexual reproduction. It 
is remarkable that a sexually reproducing population is only able to adapt faster than an 
asexually reproducing population if the recombination factor r is significantly larger than 
1. Further more we note that the threshold value of rthr is slowly increasing with increasing 
mutation rate. This is, of course, to be expected since the asexual population is able to 
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adapt faster with increasing mutation rate. 

It is clearly difficult to relate in a detailed way the mutation rate Pmut in this model 
to real biological mutation rates. Fig. |H1 shows that this is not an essential problem for 
the interpretation of our results since the same qualitative behaviour is found for all values 
of the mutation rate Pmut- Namely, recombination factors r significantly in excess of 1 is 
needed for all values of the mutation rate to make the sexual population adapt faster than 
the asexual population. 

The results in Fig. El can be understood in terms of how the velocity of the centre of 
mass of the population depends on the recombination factor for the two types of reproduc- 
tion. This is illustrated in Fig. |31 By construction the velocity of the asexually reproducing 
population is independent of the recombination factor. The velocity of a sexually reproduc- 
ing population increases rapidly with increasing recombination factor, but decreases with 
increasing maximum distance dmax between two reproducing parents. This finding may 
be understood in the following terms. Consider a reproduction event involving a most ffi 
individual, i.e. an individual on one of the positions, say S*, in genotype space with the 
highest value of the ^x-coordinate among the genotypical positions occupied at the present 
instant in time. I.e., the positions in genotype space with 5*^ > Sx are all unoccupied. An 
increase in the recombination factor leads to an increase in the range of possible S^^-values 
of the offspring's genotype and thereby assists the displacement of the population to large 
values of S^. Consider again this most fit individual S*. The average reproduction rate of 
S* is obtained from Eq. Q by averaging over all potential individuals present within the 
distance dmax- The larger the value of d^ax the smaller the Sx values of the possible mates. 
Small values of lead to small reproduction values according to Eqs. (|T)), (0) and (jH). 

The results discussed in this section are all for the case where no ffiness holes are 
present. In this case the simulations can readily be compared with the iterative master 
equation for the population density in genotype space n(S,t). These are cumbersome to 
write down but can easily be solve by numerical methods. In this way we have confirmed 
the ensemble averaged results described above. 
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4 The effects of holes 



In reality many totally unviable genotypes exist. We include this in the model by introducing 
a random selection of genotype positions which are unable to reproduce. Individuals can 
arrive at one of these positions as a result of mutation from a nearby site. Individuals on 
such a site, say Sh, may be selected for reproduction attempts. In both the sexual and 
asexual case this reproduction attempt can, however, not lead to offspring since (t>{Sh) = 0. 
The individuals on the trap sites leads to an overall fall in the reproduction rate on other 
sites through their contribute to the —nN{t) term in W, see Eq. |2l This is only natural as all 
individuals existing at a given moment in time represent a demand on the carrying capacity 
of the system. One might expect that the sexually reproducing population would be more 
likely to find a path through this rugged fitness landscape than an asexually reproducing 
population. After all, the sexually reproducing population may jump over holes in the fitness 
landscape when dmax > 2. In fact, as soon as dmax > the sexual population can pass to 
next nearest neighbour sites of already occupied positions in genotype space, whereas the 
asexual offspring can at most move to a nearest neighbour site when a mutation occur. This 
enables the sexual population to follow paths through the rugged fitness landscape which 
are inaccessible to the asexual population. 

In Fig. Elwe show as function of the density of holes the probability that a population 
is able to find an adaptive path through the rough fitness landscape in the course of a fixed 
number of time steps T = 10^ for mutation rate Pmut = 0.01. We present results for two 
cases: A) dmax = with r = 1.270 and B) dmax = with r = 1.408. In both cases r 
is calibrated such that, in the absence of holes, the sexual population moves with the same 
velocity through genotype space, as the asexual population with pmut = 0.01 does. 

The simulations clearly show that, once holes are present, the asexual reproduction 
in general is better suited to move through the rough fitness landscape, and that large values 
of dmax, rather than allowing the sexual population to more easily find a way around the 
traps, has the opposite effect. This behaviour is explained by the fact that the sexually 
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reproducing population more often place offspring in the holes than the asexual population 
does. 

5 Discussion 

It is interesting to relate our results to the width of the population along the x-axis in 
genotype space and to the width of the distribution of the reproduction probabilities. We 
find that the population that is able to adapt most rapidly always has the broadest distri- 
bution of reproduction probabilities Poff- This is entirely consistent with the Weismann's 
idea that greater variation enables adaptation to occur faster because natural selection has 
more choice to act upon [Burt, 2000| . What is not consistent with Weismann's hypothe- 
sis is that sexual reproduction does not always lead to the larger variation. In Fig. IHlwe 
show the average velocity of adaptation of the sexual and asexual population as function of 
the standard variation 0"^^^.^. of the reproduction rate of the respective population for many 
different combinations of control parameters. 

We notice that adaptive velocity and cr^^^^ are linearly related and that only for very 
substantial recombination factors is the sexually reproducing population able to adapt faster 
than the asexual. When this happens the sexual population does have the greatest variation 
in Poff- Recombination factors r much larger than one seems of little biological relevance 
since r > 1 frequently leads to offspring which may be reproductively isolated from one of 
its parents in the sense that the distance to one of the parents may exceed dmax- Thus, in 
the model considered here a sexual population is unable to adapt faster than an equivalent 
asexual population except for, what seems to be, biologically extreme choices of parameters. 

The reason we find that Weismann's hypothesis is unable to explain the superiority 
of sexual reproduction is not that the idea that greater variation leads to faster adaptation 
is wrong. The problem is that sexual reproduction events involving the most fit individuals 
will typically involve less fit mates. This hinders the individuals with the best adapted 
genotype to spearhead the adaptation of the sexual population. 

Of course the Weismann hypothesis is not the only reason why sexual reproduction 
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might be evolutionary superior to asexual reproduction. Kondrashov discuss a large number 
of different hypotheses and emphasis in particular the ability of recombination to eliminate 
deleterious mutations from the genome [Kondrashov, 1993| . Perhaps no single general mech- 
anism can be identified as the most important for the evolution and maintenance of sex (see 
e.g. [Hurst and Peck, 1996| ) in which case the explanation may have to rely on case specific 
mechanisms (see e.g. [Peck and Waxman, 2000l [Doncaster et al., 2000| ). 

We have studied Weismann's hypothesis in it simplest form and have been forced 
to conclude that within the framework of individual genotypical fitnesses combined, in the 
sexual case, to give the fitness of the reproducing pair, the extra variation in genotype 
space of the sexual population is in realistic cases not able to produce the needed excess 
in reproductive fitness to make the sexual population able to out compete the asexual 
population. The sexual population's greater variation of genotypes tends to be neutralized 
by the necessity of the most fit individual to mate with individuals of lower fitness. 

We are greatful to A. Burt and C. Godfray for comments. We thank K. Dahlstedt 
and G. Pruessner for assistance with prrgramming and P. Anderson for lingusitic advice. 
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Figure 1: Ellipse with r = lA. Si and 5*2 denote the sites of the two parents. The distance 
a is the major axis defined in the text. 
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Figure 2: The peak in the population density moving to the right along the axis in 
genotype space. The data is for a sexually reproducing population with pmut — 0.01, dmax — 
V2 and r = 1.27. The time between each small peak is equal to 6000 time steps. 
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Figure 3: The competition between sexual and asexual reproduction for two different values 

of dmax- The signatures are: A corresponds to d^ax = and □ to d^ax = V^- For a 
given dmax, sexual reproduction most efficiently response to the selective pressure above the 
relevant curve, below the curve asexual reproduction is superior. 
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Figure 4: The adaptive velocity of a sexually reproducing population for different values of 
the recombination factor and different values of d^ax = and \/5, □, respectively. The 

straight fine represent the velocity (229x 10~^ square/time step) of an asexually reproducing 
population for the same value of the mutation rate Pmut = 0.01. The two intersection points 
mark the minimum recombination factors at which sexual reproduction is superior to asexual 
ditto. 
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density of holes 

Figure 5: The fraction of the population that on average is able to pass through a rough 
fitness landscape with traps. The signatures are as follows: o asexual population, A sexual 
population with dmax = and r = 1,27, □ sexual population with dmax = and 
r = 1.41. In all cases Pmut — 0.01 and each data is averaged over 500 reahsations. 
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Figure 6: Average adaptation velocity as function of the standard deviation of the distribu- 
tion of reproduction rates. The data are for p^^t — 0.01 and 0.03, dmax — and y/5 for a 
range of different recombination factors r. 
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